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The use of spikes to carry information between brain areas implies complete 
or partial synchronization of the neurons involved. The degree of synchroniza- 
tion reached by two coupled systems and the energy cost of maintaining their 

synchronized behaviour is highly dependent on the nature of the systems. For 
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■T^lj- . non-identical systems the maintenance of a synchronized regime is energetically 

QQ ' a costly process. In this work, we study conditions under which two non-identical 



o 

(N 



electrically coupled neurons can reach an efficient regime of synchronization at 

f~^ , low energy cost. We show that the energy consumption required to keep the 

synchronized regime can be spontaneously reduced if the receiving neuron has 

adaptive mechanisms able to bring its biological parameters closer in value to 
- 1 — I 
^^ ' the corresponding ones in the sending neuron. 
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1. Introduction 

Neural information processing depends on communication between neurons 
that involves partial or complete synchronization in their signalling. When a 
neuron receives signals from its neighbours its firing regime undergoes changes 
that lead to a modification of its information capacity as well as its average en- 
ergy consumption. The degree of synchronization reached by two neurons con- 
ditions their capacity to exchange information and the energy cost to keep their 
synchronous signalling activity. Two structurally similar coupled systems can 
reach different degrees of synchronization depending on whether they are identi- 
cal or not. Identical systems reach complete synchronization spontaneously be- 
yond a given value of the coupling strength and the energy required to maintain 
this completely synchronized regime is zero. However, between non-identical 
systems complete synchronization never occurs spontaneously without cost and 
an average nonzero flow of energy is required to maintain a completely synchro- 



nized regime 
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In this work we investigate, in terms of energy consumption, 
how an efficient degree of synchronization between two non-identical electri- 
cally coupled Hindmarsh-Rose neurons can be reached. We show that feedback 
synchronization at sufficiently large coupling force between two non-identical 
neurons creates appropriate conditions for efficient actuation of adaptive laws 
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able to make the neurons approach each other their biological parameters values 
in order to decrease the energy consumption. In the next section, we describe 
the four dimensional Hindmarsh-Rose neuron model used to represent the dy- 
namics of real neurons, and we report an energy function associated to that 
model. The analysis of energy consumption during the synchronization process 
between the two neurons, and the description of the adaptive laws are performed 
in Section 3. Numerical results are showed and analyzed in section 4. Finally, 
in Section 5 we collect the main conclusions of this work. 

2. The Hindmarsh-Rose neuron energy 

In this paper, we represent a single neuron by the four-dimensional Hindmarsh- 
Rose neuron model described by the following equations of movement [2|, y, y , [Sj : 

i = ay + bx^ — cx'^ — dz + ^/, 
y^e-fx'^-y- gw, 
z ~ m{-~z + s{x + /i)), 
w = n{—kw + r{y + I)), 

where a, b, c, d, S,, I, e, /, g, m, s, h, n, fc, r, and I are the parameters that govern 
the dynamics of the neural system. The variable a; is a voltage associated to the 
membrane potential, variable y although in principle associated to a recovery 
current of fast ions has dimensions of voltage, variable z is a slow adaptation 
current associated to slow ions, and variable w represents an even slower process 
than variable z 'Gj. / is a external current input. 

In a previous work, we have assigned an energy function to the Hindmarsh- 
Rose model given by Eq. ([Ij. The procedure followed to find this energy has 
been reported in detail in [7|. This energy function H{x.) is given by 

H^^-ilfx^ + ^^^^^^x^ + ay') ^^^ 

+ ^{—^{msd — gnrjz'^ — 2dyz + 2gxw) 

In the model time is dimensionless and every adding term in Eq. ([2]) has dimen- 
sions of square voltage, so function H is dimensionally consistent with a physical 



energy as long as parameter p has dimensions of conductance. In this paper we 
fix parameter p to the arbitrary value p = —IS. We have adopted a negative 
sign for this parameter so that the outcome of the model will be consistent with 
the usual assumption of a demand of energy associated with the repolarization 
period of the membrane potential and also with its refractory period. 
The energy derivative H is given by 
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which is also dimensionally consistent with a dissipation of energy. 

The energy derivative given by Eq. ([3]) represents the net energy variation 
over time of the energy of the neuron. For an isolated neuron, whose dynamics 
is confined to an attractive region of the state space, its average net variation 
of energy is zero what indicates that the energy the neuron obtains through 
the membrane is perfectly balanced by its dissipation of energy. Nevertheless, 
when the neuron is forced to synchronize to another neuron its balance of energy 
through the membrane is broken. In this case, the coupling device must coun- 
terbalance the net energy flow through the membrane in order to lead to a zero 
global average variation. For the particular system given by Eq. ([5]), the net 
variation of energy at the coupling device is given by the following expression 
according to Ref. ;8| , where the first vector stands for the gradient of the energy 
function given by Eq. ([2]) and the second term represents the coupling device 
(see Eq. (lU). T denotes transpose of a matrix. 
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For the isolated neuron of Eq. ([I}, we have computed both energy and 



energy derivative corresponding to a series of action potentials. The results are 
depicted in Figure 1. Fig. 1 (a) shows a series of action potentials (variable 
X in the model neuron). Fig. 1 (b) and Fig. 1 (c) show both energy and 
energy derivative corresponding to that action potentials. Fig. 1 (d) shows 
two action potentials. Fig. 1 (e) and Fig. [2] (f) show detail of energy and 
energy derivative associated to a train of two action potentials. The energy 
values are negative but have been depicted in Fig. 1 (and later in figures 2 and 
5) with a positive sign for a better appreciation of the energy demand when 
generating a spike. For each action potential it can be appreciated that the 
energy derivative is first negative, dissipation of energy while the membrane 
potential depolarizes during the rising phase of the spike, and then positive, 
contribution of energy to repolarize the membrane potential during its falling 
phase. During the refractory period between the two spikes the energy derivative 
remains slightly positive, still demanding energy, until the onset of the following 
action potential. 
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Figure 1: (a) Action potentials, (b) energy and (c) energy derivative for the Hindmarsh-Rose 
model neuron, (d), (e) and (f) Details of the action potential, energy and energy derivative 
associated to two spikes 



In the next section, we analyze the average of energy and energy derivative 
of two nonidcntical coupled Hindmarsh-Rose neurons. We consider that, in a 



first stage, the coupled neurons have been forced to synchronize with each other, 
then we initiate an adaptive process that adapts the mismatched parameters of 
the receiving neuron to the one of the sending neuron. 

3. Adaptive mechanism 

Gap junction channels permit the intracellular electrical potentials of two 
neurons to directly connect together and are usually referred to as electrical 
synapses. They are frequent when two or more neurons are coupled together 
and play an important role in the synchronization of cellular events. In par- 
ticular, they are efficient in transmitting information and in synchronizing the 
information of groups of neurons [9| . Although symmetry is an expected prop- 
erty for an electrical synapse, asymmetric gap junctions have also been reported 
I I 

in the literature [10]. In this section we analyze the energy requirements in- 
volved in achieving a complete synchronization between two initially noniden- 
tical Hindmarsh-Rose neurons. We first consider an asymmetric unidirectional 
coupling scheme in which only the first variable of the receiving neuron is af- 
fected by the coupling according to the following system of equations: 



(5) 



it = ayi + bxf - dxf - dz, + ^li + k,{xj - x^), 
yi = ei- fixf ~yi~ gwt, 

Zi = m(-Zi + s{xi -\- h)), 

Wt = n{-kwi + r{yi + /)), 
where ki = and k2 > is the coupling strength. Note that the coupling 
affects only the membrane voltages X2 of the receiving neuron, i, j — l,2;i ^ j 
are the indexes for the neurons. 

We suppose that the sending neuron is set in the chaotic spiking-bursting 
regime corresponding to a constant external current /i = 3.024, with the stan- 
dard parameter values a = 1, b = 3.0{mV)~'^, c\ — \[mV)^'^ , d = Q.99Mil, 
^ ^ IMn, ei = l.OlmV", /i = 5.0128{mV)-\ g = 0.0278Am, m = 0.00215, 
s = 3.966^^5", h = 1.605my, n = 0.0009, k = 0.9573, r = 3.0^^5*, I = 1.619mT/. 



The receiving neuron is initially set to its quiescent state at a low value 
I2 = 0.85 of its external current, the parameters C2, 62 and /2 are mismatched 
as follow: C2 = 0.95(mF)-2, £2 = 0.85mT/, /2 = 5.1128(ml/)-i. All others 
parameters are the same for both neurons. 

Under these conditions the coupled neurons are not identical and there- 
fore synchronization does not spontaneously occur at any value of the coupling 
strength but, rather, it must be strongly enforced through the establishment of 
large values of the coupling strength which will involve high energy consump- 
tion. The energy consumption required to maintain the synchronized regime 
can be spontaneously reduced if the receiving neuron is flexible enough as to 
adapt its parameters through an adequate adaptive law in order to reach the 
nominal value of the sending neuron parameters. Ideally, if the neurons become 
identical their joint dynamics is attracted toward a regime of zero error in the 
variables. This asymptotical limit regime of identical synchronization occurs 
with zero net average exchange of energy with the environment. 

In a previous work 'Sj we deduced adaptive laws that permit to any family of 
homochaotic coupled systems to adapt their structures in order to approach each 
other provided they are previously forced to a certain degree of synchronization. 

The coupling scheme given by Eq. ([5]) can be written, in general, as follow: 

xi=/(xi,p), 

(6) 
X2 = /(X2,q) +K(xi -X2), 

where xi,X2 G 5ft" indicate the states of the coupled neurons, p and q stand 

for the parameters of the sending and receiving neurons, and that K e K" is 

a diagonal matrix representing the coupling strength with entries k; > 0,i = 

1, ...,n. In the particular case where only the first variable of the vector X2 is 

affected by the coupling, we have ki = k and k^ = 0, i = 2, ..., n. 

If the coupling strength k is large enough as to make the errors in the 

variables C = X2 — Xi small, an operational law that adapts the parameters of 



the receiving neuron to the ones of the sending neuron is given by [8| 
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c = 



I V ^. /(xi,p) J /^-j 

where C'' ^ Q ^ P denotes the vector of parameter errors, and the summation is 
over every component of the vector field /. The above law is general and can be 
used to find specific adaptive laws to any kind of homochaotic systems provided 
they are coupled through a feedback scheme of large enough coupling strength. 
Based on the above adaptive mechanism, and according to Eq. ([T]), the 
adaptive laws that govern the dynamics of the mismatched parameters of the 
receiving neuron of the system of Eq. ([S|), are given by the following equations: 



h = -i{x2 - Xi) 

C2 = X\{X2 - Xi) + 3xl{x2 - Xi)^ 

62 = -(2/2 ~y\) 

h = x\{y2 - yi) + 2a;i(a;2 - xi)[y2 - y\) 



(8) 



4. Numerical results 

Numerical results have been performed simulating the system of Eq. (|5|) 
over 50000 time units. 

In a first phase, we force the receiving neuron to replicate the behavior of 
the sending neuron through a progressive increase of the coupling force k to a 
value large enough to ensure a regime close to complete synchronization. When 
k reaches the appropriate value we keep it fixed for a second phase in which 
we begin the process of parameters adaptation. During the forced synchroniza- 
tion the parameter k has been increased linearly from to 30, value that has 
remained fixed during the adaptation phase. The registrated values correspond- 
ing to energy and energy derivative have been averaged over a convenient length 
of time in order to avoid large fluctuations. The dissipated energy has been av- 
eraged over five units of time, while the proper energy has been averaged over 
ten units of time. The energy at the coupling device is calculated according to 
Eq. Q. 
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Figure 2: In (a) Average energy over five units of time of sending and receiving neurons. In 
(b) average per unit time of the energy variation over ten units of time. Adaptation begins 
at t = 30, 000. 



Figures [2][a) reports the average values of energy per unit time of sending 
and receiving neurons. We can see that in a first stage {t < 30,000), the 
receiving neuron is forced to synchronize with the sending neuron, and shows a 
decreasing pattern of average energy due to the gradual increase of the coupling 
force. In this regime, the receiving neuron is continuously demanding energy 
from the coupling device, leading to an increase of the energy variation through 
the membrane as we can see in part (b) of the Fig|21 The energy demanded by 
the receiving neuron must be supplied by the coupling device. After adaptation 
takes place, the two neurons become structurally close each other, and enter in a 
completely synchronized regime of zero error in the variables as shown in Fig. 21 
This regime corresponds to a balanced exchange of energy between the neuron 
and its environment corresponding to a zero value of the energy derivative. 

The energy consumption of the sending neuron remains constant because its 
dynamics is not affected by the coupling, this energy is perfectly balanced by 
the income of energy it receives through the membrane, i.e. its energy variation 
is Hi = 0. 

In Fig. [3] we report the results of the adaptation mechanism that has been 
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Figure 3: In (a)-(d) adaptation results of the mismatched parameters of the receiving neuron. 

implemented according to Eq. ([S]) to adjust the mismatched parameters of the 
receiving neuron to those of the sending neuron. As it can be appreciated, 
a correct adaptation of the accessible parameters is achieved. In Fig. |4l we 
can see that during the first phase, before the adaptation process begins, the 
two neurons have already experienced a notably decrease in the synchronization 
error which correspond to a high level of synchronization. This synchronized 
signalling implies that information is efficient in the sense that it is being trans- 
mitted between both neurons at high signal to noise ratio. However, as Fig. 
m^b) shows, to keep this synchronised regime implies a high energy cost to 
the coupling device. In other words, maintaining the synchronized regime for 
communication has an additional energy cost over the basic energy cost of the 
normal neuron's metabolism. This additional cost is direct consequence of the 
two neurons being non-identical. Figure [^Jb) shows in a second phase that, if 
the receiving neuron adapts its biological parameter values, the two neurons 
synchronize with zero net flow of energy from the coupling device. That is, the 
efficiency of the information transmission is slightly increased and transmission 
itself has no energy cost. 

As bidirectional couplings between neurons are frequent, we have also con- 
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Figure 4: Synchronization errors in state space variables. 

sidered a bidirectional electrical coupling corresponding to ki = k2 = k in Eq. 
([S]). In this case, the mutual interaction induced by the coupling affects both 
the variable xi and X2 of the sending and receiving neurons. As in the previ- 
ous experiment, we analyze the energy consumption required to maintain the 
synchronized regime between the coupled neurons, and the energy average that 
they receive from the external source, i.e., the coupling device. 

We suppose that the dynamics of the sending neuron is only affected by 
the coupling but not by the mechanism of parameters adaptation, so that all 
its parameter values are kept fixed whereas the mismatched parameters of the 
receiving neuron are governed by the adaptive laws given by Eq. ([5]) . 

Figure [5] shows the average values of energy and energy derivative per unit 
time of the sending and receiving neurons. 

As it can appreciated, before adaptation starts (f < 30,000), and as soon 
as the coupling device is connected, the coupled neurons begin to adjust their 
dynamics to achieve a common behaviour. At this stage, the sending neuron, 
as shown in Fig. [SJb) , undergoes an increase in its average energy which corre- 
sponds to an increase of its dissipation of energy (negative values of its energy 
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derivative). On the other hand, the receiving neuron displays an opposite pat- 
tern corresponding to a decrease in its average energy (see Fig. [S]) and an 
increase in its demand of energy to be supplied by the coupling device. After 
adaptation takes place, the neurons reduce notably their average energy values 
to reach a value of about 50 corresponding to the average energy of an isolated 
neuron (see Fig. Elja)). At this stage the energy variation of both neurons is 
zero. 




Figure 5: In (a) energy average over five units of time of the sending and the receiving neurons. 
In (b) average per unit time of the energy variation over ten units of time of the sending and 
the receiving neurons. Adaptation begins at t = 30, 000. 



5. Conclusion 

Signalling in synchrony is a normal way to propagate information between 
neurons. However, for non-identical neurons the cost of signalling in synchrony 
is superior to the mere metabolic energy cost that the neurons would have if 
they were signalling independently. Maintaining a synchronized regime between 
two non-identical neurons requires energy that must be provided by the coupling 
biological device that forces that synchrony. In contrast with non-identical neu- 
rons, identical neurons could reach a spontaneous synchronized regime at zero 
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maintenance energy cost. 

This paper analyzes the energy cost of maintaining a synchronized regime 
between two electrically coupled Hindmarsh-Rose neurons. We have consid- 
ered a sending neuron always signalling in a chaotic regime and a non-identical 
receiving neuron, with some mismatch in its parameters, initially set in its qui- 
escent state. Under these conditions, to keep the synchronized regime requires a 
net flow of energy that can be costly to maintain. We show how this energy cost 
reduces notably when the receiving neuron is able to adapt its structure in order 
to approach the dynamics of the sending neuron. Biological structures are par- 
ticularly flexible in adapting their parameters and the mechanism of adaptation 
introduced in this work could make some of the required collective behaviours 
between groups of neurons energetically less costly. 
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